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Abstract 

The formation of the adatom adsorption structure in dynamic force microscopy experiment is shown as 
a result of the spontaneous appearance of shear strain caused by external supercritical heating. This 
transition is described by the Kelvin- Voigt and Maxwell-type equations for a viscoelastic medium and by 
the relaxation equation for temperature. These equations coincide formally with the Lorenz system. The 
steady-state values of the shear strain and stress, and the sample temperature are found. The critical 
temperature of the tip increases with growth of the effective value of the sample shear modulus and 
decreases with growth of its typical value linearly. 
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1. Introduction 

Nowadays due to large scientific and practical importance the phenomena, taking place on the sample 
surface at interaction with tip of dynamic force microscope, e.g., atomic force microscope (AFM) and 
friction force microscope, attract more and more attention (see the reviews in and the 

literature cited therein). Particularly, the experimental and theoretical data are obtained on structural 
instabilities, phase transformations, plastic dislocation, neck and adatoms structures formation 0, a B 
a HqI, U, 12, 13||. These processes are characterized by hysteresis of dependencies of adhesion force and 
potential energy surface on the tip-surface distance [10|, [14|, |l5| and by hysteresis of the sample stress vs 
strain curve 
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Since the nature of such phenomena remains poo rly understood the basic goal of the present study is 
the construction of a qualitative nonlinear model [l7|, [l8|, [l9|, |20|, |2l|, |22] describing the hysteresis processes 
which occur on the germanium surface during interaction with the AFM tip 10J. Here the macroscopic 
continuum mechanics models |23|, |24| are supposed to be still applicable to the atomic length-scales, 
where discrete atomistic interactions become significant 0, [25|, 26]. However, a total explanation of 
studied macroscopic phenomena requires the consideration of the microscopic processes. In the presented 
approach the formation conditions of the adatom adsorption structure are defined on the semiconductor 
surface due to both thermal and deformation effects. The total set of freedom degrees is considered as 
equivalent variables. The adatom configuration formation is described analytically as a result of the 
self-organization caused by the positive feedback of shear strain and temperature on shear stress on the 
one hand, as well as the negative feedback of shear strain and stress on temperature on the other hand. 
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This study is based on the assumption that stress relaxation time diverges because the shear modulus 
vanishes at the point of transition. 

The paper is organized as follows. In Section [2] the self-consistent Lorenz system of the governed 
equations is written for a viscoelastic semiconductor approximation subject to the heat conductivity. 
The adatom structure formation is shown in Section [3] to be supercritical in character (has the second 
order) when the effective shear modulus of the germanium does not depend on the strain value; it then 
transforms to a subcritical mode with this dependence appearance (Section H]). In these sections the 
steady-state values of shear strain and stress, as well as temperature are determined also within adiabatic 
approximation. Using such limit, the synergetic potential is obtained, that is the analog of thermodynamic 
potential, from basic evolution equations. Section [5] contains a short conclusions. 



2. Basic equations 



Let us start from the supposition that the relaxation behavior of the shear component e of the strain 



tensor in a viscoelastic semiconductor is governed by the Kelvin- Voigt equation [23 



e = -e/r £ + a/r) E 



(1) 



where r E is the Debye relaxation time and rj E is the effective shear viscosity coefficient. The second term 
on the right-hand side describes the flow of a viscous liquid caused by the corresponding shear component 
of the stress a. In the steady state, e = 0, we obtain the Hooke-type expression a = G E e, G E = t] e /t e . 

The next assumption of our approach is that the relaxation equation of the sample shear stress a has 
a similar form to the Landau-Khalatnikov equation 



19. 27. 21 



r a a = —a + Ge. 



(2) 



Here the first term on the right-hand side describes the relaxation during time r a = t]/G determined by 
values of the shear viscosity 77 and modulus G. In the stationary case & = the kinetic equation fl2]) is 
transformed into the Hooke's law 

a = Ge. (3) 

Substituting de/dt instead of e/r a in Eq. reduces it to a Maxwell- type equation for a viscoelastic 
matter |24( . Eqs. ([I]) and ([2]) give a new rheological model because they can be transformed to a differential 
equation of second order with respect to strain or stress. 

Note that effective values of viscosity r] e = t e G £ and modulus G E = t] e /t e do not coincide with the real 
values 77 and G. Physically such difference is conditioned by that the Maxwell-type equation fl2]) is not 



equivalent to the Kelvin- Voigt equation ([[]) [23J, |24|, |2l|. As is known the values G E , 77, r\ E depend very 
weakly on the s ample temperature T, while the real shear modulus G vanishes, when the temperature 

Further, the simplest approximate temperature dependencies are used: 

G (T/T C -1), (4) 



decreases to T r [28L 29, 30, 31 



G E (T), r)(T),r] E (T) = const, 



G 



where G = G(T = 2T C ) is the typical value of modulus. 



According to the synergetic concept [17|, |20|, l21|, |3(J, |31|, |32j, [33J for completing the equations system 
([T]) and ([2]), which contains the order parameter e, the conjugate field a, and the control parameter 
T, we should deduce a kinetic equation for the temperature. This equation can be obtained using the 
basic relationships of elasticity theory stated in § 31 in Ref. j24j]. Thus, it is necessary to start from the 
continuity equation for the heat 5Q = T5S: 



TS = -Vq. 



(5) 
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Here the heat current is given by the Onsager equation 



q = -kVT, (6) 
where k is the heat conductivity. In the elementary case of the thermoelastic stress the entropy 

S = S (T) + Kae° (7) 
consists of the purely thermodynamic component So and the dilatation: 

e° = e°l e° = a(T- T ) , (8) 

where a is the thermal expansion coefficient, To is the equilibrium temperature, I is the unit tensor 
and K is the compression modulus (see § 6 in (24|). In the considered situation we should transfer 
from the dilatational component Kae° to the elastic energy —as/T of the shear component divided by 
temperature (here the minus sign takes into account the connection T5S = p5V =>■ — a5e at Sq = const, 
which is caused by the contrary choice of the pressure p and the stress a signs). As a result, Eq. (jSJ) 
possesses the form 

TS (T) -ae= kV 2 T. (9) 

Taking into account the equality (k// 2 )(T 6 — T) « kV 2 T (I is the scale of heat conductivity, T e is the 
AFM tip temperature) and the definition of heat capacity c p =TdSo/dT, Eq. ([9]) assumes the form: 



c p T 



K 



I 



2 {T e -T) + ae. 



(10) 



Substituting here expression for the e from Eq. (TTJ we obtain the term <7 2 /r] £ . It describes the dissipative 
heating of a viscous liquid, flowing under influence of the stress a, that can be neglected in the case under 
consideration. On the other hand, the process of an AFM tip moving into contact with the surface has 
the following peculiarity. It is necessary to consider the thermal influence of the tip whose value T e is not 
reduced to the Onsager component and is fixed by external conditions. In view of these circumstances 
the square contribution of the stress is supposed to be included in T e . The obvious account of this term 
leads to a significant complication of the subsequent analysis, though it results in a renormalization of 
quantities. Therefore component T e in Eq. (flQl) is assumed to be constant for our further consideration. 
It is convenient to introduce the following measure units: 



a s = (c p r] £ T c /T T ) 



1/2 



e s = a s /G £ , 



for the variables a, e, T, respectively (tt = l 2 c p /n is the time of heat conductivity). Then, the basic 
equations (PQ), fl2}, and ffTOj) take the form: 



r £ e = — e + a, 

r a & = -a + g(T - l)e, 

r T f = (T e — T) — CTE, 



where the constant 
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is introduced. Equations (|T2|) - 

to describe the thermodynamic phase and the kinetic transitions [20|, [21 



have a similar form that of the Lorenz scheme QJ 

3c 




(12) 
(13) 
(14) 



(15) 



which allows us 
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3. Continuous transition 



In general the system f|T2|) - ( Tl4j) can not be solved analytically, therefore we use the following adiabatic 
approximation: 

Ta < r £ , r T < r 6 . (16) 

This implies that in the course of matter evolution the stress a(t) and the temperature T{t) follow the 
variation of the strain e(t). The first of these inequalities is fulfilled because it contains the macroscopic 
time t £ and the microscopic Debye time r CT a/c~10~ 12 s, where a ~ 1 nm is the lattice constant or the 
intermolecular distance and c ~ 10 3 m/s is the sound velocity. The second condition ([16]) can be reduced 
to the form 

/ < L, (17) 
where the maximal value of the characteristic length of the heat conductivity 
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the thermometric conductivity % = k/c p , the effective kinematic viscosity v £ = i] e /p and the sound 
velocity c £ = (G e /p) 1//2 are introduced (p is the medium density). Then, we can put the left-hand sides 
of Eqs. (|T3|) and (|T4|) to be equal to zero. As a result, the stress a and the temperature T are expressed 
in terms of the strain e: 



1 + ge* ' 
T e - 1 

T = 1 + -i— j. 20 

1 + ge l 



In accordance with Eq. ()2Up . in the important range of values of the parameter T e > 1, the temperature T 
decreases monotonically with increasing strain e from the value T e at e = to (T e +l)/2 at e = e m = sfvfg- 
Obviously, this decrease is caused by the negative feedback of the stress and the strain on the temperature 
in Eq. (1141) . that is explained by the Le Chatelier principle for this problem. Really the reason for the 
formation of adatoms adsorption structure is the positive feedback of the strain and the temperature on 
the stress in Eq. ( 1131) . Hence the increase in the temperature should intensify the self-organization effect. 
However according to Eq. (I14p the system behaves so that the consequence of transition, i.e., growth 
of the strain, leads to a decrease in its reason (temperature). Equation (fT9|) . expressing the stress in 
terms of the strain, has the linear form of the Hooke's law at e C e m with the effective shear modulus 
G e f = g (T e — 1). At e = e m the function a(e) has a maximum and at e > e m it decreases that has no 
physical meaning. Thus, the constant e m = y/l/g gives the maximal strain. The increase in the typical 
value of the modulus Go leads to a decrease in the maximal strain e m and an increase in the effective 
modulus G e f whose value is proportional to the characteristic temperature T e . 



Substituting Eq. (HM into Eq. f[T2"j) we obtain the Landau-Khalatnikov-type equation [19|, [27 

r £ e = -dV/de, (21) 

where the synergetic potential has the form 

V = ±[e 2 + (l-T e )ln(l + ge 2 )]. (22) 



At steady state the condition e = is fulfilled and the potential ( 1221) acquires a minimum. When the 
temperature T e becomes smaller than the critical value 

T d} = l + g- 1 ; g = G /G £ <l, G £ = V£ /t £ , (23) 



4 



this minimum corresponds to e — 0, i.e., the adatom adsorption structure is not realized. In the reverse 
case T e > T c q, the stationary shear strain has the nonzero value 

e =(T e -{l+g- l )) l/2 (24) 

which increases with T e growth in accordance with the root law. This causes the formation of the adatom 
configuration. Equations ( fl~9l) and ( J201) give the stationary values of stress and temperature: 

a = e , T = l + g-\ (25) 

Note that, on the one hand, the steady temperature To coincides with the critical value (|23|) and, on 
the other hand, its value differs from the temperature T e . Since T c q is the minimal temperature at 
which the formation of the adatom adsorption structures can be observed, the above implies that the 
negative feedback of the stress o and the strain e on the temperature T (see last term on the right-hand 
side of Eq. ffT4j) ) decreases the sample temperature so much that only in the limit does it ensure the 
self-organization process. At steady state the value of the shear modulus is 

G s = G £ . (26) 

The two cases can be marked out by the parameter g = Gq/G £ . In the situation g ^> 1, meeting the 
large value of the modulus Go, Eqs. ([23]) - f )25|) take the form 

e = (T e - T Q = T c0 = 1. (27) 

This corresponds to the "solid (fragile)" limit. The opposite case g <C 1 (small modulus Go) meets the 
"strongly viscous liquid" 

6 = (T e - 9' 1 ) 1 ' 2 , To = T c0 = g- 1 = G £ /G . (28) 



4. Effect of deformational defect of modulus 

The Kelvin- Voigt equation ([T]) assumes the use of the idealized Genki model. For the dependence 
a(e) of the stress on the strain, this model is described by the Hooke's expression a = G e e at e < e m 
and the constant a m = G e e m at e > e rn (a m , e m are the maximal stress and strain, a > a m results in 
viscous flow with the deformation rate e = (a — a m ) /r] £ ). Actually, the <r vs e dependence curve has 
two regions: first one, Hookean, has a large slope corresponding to the shear modulus G £ , followed by 
the more gently sloping section of the plastic deformation whose tilt is defined by the hardening factor 
< G e . Obviously, such a picture means that the shear modulus, introduced in Eq. (DQ), depends on the 
strain value. Let us use the simplest approximation 

Gs(e) = + 1 ^ e "% , (29) 
1 + (e/e p ) 

which describes the above mentioned transition of the elastic deformation mode to the plastic one. It 
takes place at a characteristic value of the strain e p , which is smaller than e s (otherwise plastic mode is 
not realized). Note that an expression of the type Eq. ff29l) was proposed, for the first time, by Haken 



17] describing the rigid mode of laser radiation. It is used |20|, |32|, |33[ to describe the first-order phase 



transition, and Eq. ( )29l) had contained the square of the ratio e/e p (so the V vs e dependencies in Ref. [20 



32l . l33| and Eq. ( 1301) have the even form). In description of the structural phase transitions of a liquid 



the third-order invariants, breaking the specified parity, are present [18] . Therefore in study [2JJ, |3C| [31 
in approximation (|29|) we used the linear term e/e p , instead of the square one (e/e p ) 2 . Obviously, in such 
case the V vs e dependence is not already even. 
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Within the adiabatic approximation ( TToT) the Lorenz equations (1121) - ( ED , where G e is replaced 
by dependence G £ (e), is reduced to the Landau-Khalatnikov equation ( I2TT) . The synergetic potential 
possesses the form: 



v = i £2 _^-D( i , n 

2 2 l#a 2 - # 



1 + ^ 2 



1 + 9{e/a) 

21 



+ w- g -'^) [rl ' n|rl + (£/ 



a 



- ^a" 2 ln^" 1 ^ 2 + (^/a) 2 ]]} . (30) 

Here the constant a=e p /e s < 1 and the parameter 6=Q/G e < 1, describing the ratio of the tilts for the 
deformation curve on the plastic and the Hookean sections, are introduced. At small value of temperature 
T e dependence (1301) has a monotonically increasing shape with its minimum at e = corresponding to 
the steady state of the absence of the adatom adsorption structure (curve 1 in Fig. [TJ. As shown in Fig. [T] 
at 



T°=l+6g- 1 +a 2 (l-26)+2a v /g- 1 6(l-6)(l-ga 2 ), (31) 

a plateau appears (curve 2), which for T e > T c ° is transformed into a minimum, meeting the strain e 7^ 0, 
and a maximum at e m that separates the minima corresponding to the values e = and e = Eq (curve 3) 



10j. With further increase in the temperature T e the "ordered" phase minimum, corresponding to the 
adatom adsorption configuration e = e , grows deeper, and the height of the interphase barrier decreases, 
vanishing at the critical value T c0 = 1 + g^ 1 (123]) . The steady-state values of the strain have the form 
(see Figs. □ and ED 



(^) 2 = (2^)- 1 {^(T e -l 



a 2 ) 



T y[g (T e -l-« 2 ) -6Y -Aga^O [l-g (T.-l)] j , (32) 

where the lower sign meets the stable adatom structure and the upper sign corresponds to the unstable 
one. At T e > T c0 the dependence V(e) is characteristic for the absence of the modulus defect (see curve 4 
in Fig. CD. 

It is worth noting that potential barrier inherent in the synergetic first-order transition is manifested 
only in the presence of the deformational defect of the modulus. Since latter is realized always [loj], it 
follows that a studied adatom structure formation is a synergetic first-order transition. The considered 
situation is much more complex than typical thermodynamic phase transitions. Really, in the latter 
case the steady-state value of the semiconductor temperature To is equal to the value T e fixed by a 
thermostat. In this study T is reduced to the critical value T c0 for the synergetic second-order transition 
(see section [3D . When the modulus defect is taken into account, the temperature 



1 + 7^7-^ (33) 



T e -1 

whose value is defined by a minimum position of the dependence ( l30~p . is realized. According to Eqs. ()32 
and ( |33l . the quantity T smoothly decreases from the value 

T°- 1 
T = 1 -I — — 



1 + 9 [el) 



e 



■ (2g6)- 1 [g (T c ° - 1 - a 2 ) - 6] (34) 



at T e = T c °, to 1 at T e — > oo. As shown in Fig. [31 the stationary temperature To increases linearly from 
to T c o, with T e being in the same interval and, after the jump down at T e = T c q, the magnitude Tq 
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Figure 1: Dependence of the synergetic potential on the strain at g=0.2, 8=a— 0.25 and various temperatures: (curve 1) 
T e <T c °, (curve 2) T e =T c °, (curve 3) T c °<T e <T c , and (curve 4) T e >T c0 . 




Figure 2: Dependence of the steady-state values of the strain on the temperature T e at parameters of Fig. Q] (the solid curve 
corresponds to the stable value e^, the dashed curve meets the unstable one, e m ). 
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Figure 3: Dependence of the steady-state value of the sample temperature To on the temperature T e at parameters of Fig.[T] 

smoothly decays. If the temperature T e then decreases, the stationary temperature T grows. When the 
point T c ° [Eq. ( 13T|) ] is reached, T undergoes a jump from T m [Eq. up to T c °. For T e < T c °, again the 
stationary temperature T is equal to T e . 

Since T c ° > 1, the maximal sample temperature ( l34j) is lower than the minimal temperature of the 
AFM tip (13T1) . As shown in Fig. |3] at T e > T c ° the stationary temperature To of the sample is lower than 
T 

5. Summary 

In accordance with the analysis presented above the formation of the adatom adsorption structure is 
caused by the self-organization of the shear components of the strain and the stress fields, on the one 
hand, and the sample temperature, on the other hand. Here, the strain e acts as the order parameter, the 
conjugate field is reduced to the stress <r, and the temperature T is the control parameter. The cause for 
self-organization is the positive feedback of T and e on a (see Eq. ({13]) ). According to Eqs. (]2]) and 01]), it 
is caused by the temperature dependence of the shear modulus. The formation of the adatom adsorption 
structure due to both heat and deformation effects at contact of solid surfaces can be represented as shear 
and dynamic disorder- driven melting where thermal fluctuations [iH are absent. With an allowance 
for the effective shear modulus vs strain dependence, we obtain the expressions for the temperatures 
corresponding to the absolute instability of the adatom configuration T c ° [Eq. f l31~]) ] and its stability limit 
T c o [Eq. (123]) ]. The real thermodynamic transition temperature can be determined from the equality 
V(0) = V(eo) of potentials in different phases and it is in the (T c °,T c0 ) region. According to Eq. f ]23]) 
systems predisposed to formation of adatom structure have large typical Go and small effective G e values 
of shear modulus. 

The kinetics of this transition is described by the Landau-Khalatnikov equation (|2T]) . where the 
synergetic potential has the form ( 130]) inherent in the first-order transition. In supercooled adatom 
configuration with r £ = oo, freezing of the system can occur (e — > 0) even in the nonsteady state 
OV/de ^ 0. 
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